Absorption Spectra of Astaxanthin Aggregates 
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Carotenoids in hydrated polar solvents form aggregates characterized by dramatic changes in their 
absorption spectra with respect to monomers. Here we analyze absorption spectra of aggregates 
of the carotenoid astaxanthin in hydrated dimethylsulfoxide. Depending on water content, two 
types of aggregates were produced: H-aggregates with absorption maximum around 390 nm, and 
J-aggregates with red-shifted absorption band peaking at wavelengths >550 nm. The large shifts 
with respect to absorption maximum of monomeric astaxanthin (470-495 nm depending on solvent) 
are caused by excitonic interaction between aggregated molecules. We applied molecular dynamics 
simulations to elucidate structure of astaxanthin dimer in water, and the resulting structure was 
used as a basis for calculations of absorption spectra. Absorption spectra of astaxanthin aggregates 
in hydrated dimethylsulfoxide were calculated using molecular exciton model with the resonance 
interaction energy between astaxanthin monomers constrained by semi-empirical quantum chemical 
calculation. The intramolecular vibrations of astaxanthin are modeled by a line shape function 
corresponding to two characteristic C-C and C=C stretching modes with frequencies of 1150 cm -1 
and 1520 cm -1 . The solvent is represented by a spectral density of harmonic vibrations. The 
spectral changes at increasing concentrations of water were assigned to formation of aggregates with 
decreasing exposure of the astaxanthin hydrophobic chain to the solvent. 
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I. INTRODUCTION 

Carotenoids are a widespread group of natural pig- 
ments that attracted a lot of attention during the past 
decade due to their rich photophysics (U-H)]. Their pho- 
tophysical properties are important for their functions 
in photosynthetic organisms, where they act as key 
molecules in both light-harvesting and photoprotection 
0, [H|. Based on a number of studies of monomeric 
carotenoids it is now a well-established fact that the 
strongly absorbing transition from the ground state re- 
sponsible for the characteristic colors of carotenoids is 
not their lowest energy transition. Instead, at least one 
dark singlet state, denoted as the S\ state, lies below the 
absorbing state which is denoted as S2 state, accordingly 
Moreover, recent results indicated that other dark 
states inaccessible from the ground state may be located 
within the S1-S2 energy gap [2|-|4|, 

While the excited-state properties of monomeric 
carotenoids were studied to great details by various ex- 
perimental and theoretical methods [ll-fl2| , spectroscopic 
properties of carotenoid aggregates, which are readily 
formed in hydrated polar solvents |13H18| . are much less 
understood. Two types of carotenoid aggregates can be 
distinguished according to their absorption spectra. The 
first type is associated with a large blue shift of the ab- 
sorption spectrum of a monomeric carotenoid that is also 
accompanied by a loss of vibrational structure of the 
S2 state. These aggregates are usually denoted as H- 
aggregates, in which the conjugated chains are closely 
packed oriented parallel to each other . The sec- 

ond aggregation type (J-aggregate) is characterized by a 



red-shift of the absorption spectrum, while the resolution 
of vibrational bands is mostly preserved. This aggrega- 
tion (J-type) is likely a result of a head-to-tail organi- 
zation of conjugated chains, forming a loose association 
of carotenoid molecules [TH, . Besides changes in ab- 
sorption spectra, carotenoid aggregates also exhibit sig- 
nificant changes in circular dichroism (CD) 13L excited- 
state dynamics [l4j|, and Raman spectra [20, l2l|. 

Recent experimental studies showed that carotenoid 
aggregates may play a significant role in various natu- 
ral and artificial systems. Carotenoids tend to aggre- 
gate when present in lipid bilayers, where they typically 
form H-type aggregates [22l.[23|. However, it was recently 
shown that absorption changes consistent with J-type ag- 
gregation may occur in certain carotenoproteins [2J, 1 2 ol | . 
In artificial systems, H-aggregates are often formed when 
carotenoids are deposited on surfaces. Since assemblies 
consisting of carotenoid molecules attached to conducting 
or semiconducting materials hold promise to act as pho- 
toactive species in dye-sensitized solar cells (26l - l30l |. un- 
derstanding the effects of aggregation on the structure of 
excited states is an important factor in controlling the ef- 
ficiency of such devices. Also, aggregation of carotenoids 
facilitates a homofission process [31|, in which the S\ 
state of carotenoid aggregate breaks down to form two 
triplets [2 (J . I'2ll | . When attached to an electron acceptor 
(e.g. Ti02 semiconductor), this process promises to gen- 
erate two electrons from a single absorbed photon, thus 
it is extremely attractive for the future design of devices 
aiming for conversion of solar energy to electricity. 

Carotenoid astaxanthin, which is the subject of this 
study, is known mainly as the colorant of the lobster 
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carapace where it is the only pigment of the caroteno- 
protein crustacyanin [32| . Upon binding to crustacyanin, 
absorption spectrum of astaxanthin undergoes a large red 
shift whose origin is still a matter of debate and exci- 
tonic interaction due to aggregation has been one of the 
suggested mechanisms of the red shift [33] . A few exper- 
imental studies of astaxanthin aggregation in hydrated 
solvents demonstrated that at ambient temperatures as- 
taxanthin predominantly forms H-aggregates in hydrated 
acetone or methanol [TEl [H| . J-aggregates are gen- 
erated and stabilized only upon heating up the system 
above 30 °C [l^. Yet, no theoretical studies aiming for 
explanation of aggregation-induced changes of absorption 
spectra of astaxanthin were carried out so far. 

In this work, we explain absorption spectra of astaxan- 
thin aggregates employing Frenkel exciton model for res- 
onance interaction within astaxanthin aggregates, con- 
strained by molecular dynamics simulations of an astax- 
anthin dimer. Two related models of astaxanthin aggre- 
gates are developed. For the purpose of calculating ab- 
sorption spectra of astaxanthins, we introduce a molec- 
ular exciton model. In order to estimate one of the es- 
sential parameters of the exciton model, the resonance 
coupling, and in order to constrain its values for sub- 
sequent fitting, we use a more detailed, but still sim- 
plified, semi-empirical 7r-electron molecular orbital model 
of the astaxanthin. Excitonic model with explicit treat- 
ment of fast vibrational modes of carotenoids was used 
by Spano et al. [19] to explain absorption and CD spec- 
tra of carotenoids. Here, the general approach is simi- 
lar, although we do not consider the vibrational modes 
explicitly, but describe them by energy gap correlation 
function. 

In Section [TT] we discuss the preparation of the as- 
taxanthin aggregates in different solutions. Section IIII1 
presents the lineshape method of calculation of absorp- 
tion spectrum. The intramolecular vibrations of the as- 
taxanthin are modeled by a two modes with frequencies 
u< vib = 1150 cm -1 and Lu^ b — 1520 cm -1 and the sol- 
vent is represented by a spectral density of harmonic 
vibrations (e.g. instantaneous normal modes). In Sec- 
tion IIVI we discuss the molecular exciton model, where 
the astaxanthin is represented by two electronic levels, 
(the ground state \g) and exited state |e)) standing for 
the So and the allowed S2 states of the astaxanthin, re- 
spectively. We also discuss calculation of the absorption 
spectrum of the aggregates of astaxanthins based on such 
a model. The crucial parameter, which characterizes the 
aggregates assembled after the water is added into the so- 
lution, is the resonance coupling between the members of 
the aggregate. We perform molecular dynamics simula- 
tions in Section [V] to obtain an insight into the structure 
of the astaxanthin dimer in water. The average distance 
between the molecules in this case is used to qualita- 
tively estimate the resonance coupling in Section fVl] In 
Section I VIII we compare the spectra of different aggre- 
gates of astaxanthin and fit the experimental absorption 
spectra obtained under various solution conditions. 
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Figure 1: Part A: Absorption spectra of astaxanthin in pure 
DMSO and in various DMSO/water mixtures. In all sam- 
ples, the concentration of astaxanthin is 20 (micro)M. Part 
B: Molecular structure of astaxanthin with the part of the 
7T— conjugated system used for calculations marked in light 
blue. 



II. PREPARATION OF ASTAXANTHIN 
AGGREGATES 

Absorption spectra of astaxanthin monomers and ag- 
gregates in hydrated dimethylsulfoxide (DMSO) along 
with molecular structure of astaxanthin, are shown in 
Fig. [T] DMSO was chosen due to good solubility of as- 
taxanthin in this solvent and also due to the large vis- 
cosity of DMSO which stabilizes aggregates. Also, as 
descried below, both types of aggregates can be read- 
ily formed in DMSO. Absorption spectrum of astaxan- 
thin monomer in DMSO peaks at 20280 cm" 1 (493 nm). 
The S0-S2 transition of astaxanthin monomer does not 
exhibit any vibrational structure, a phenomenon charac- 
teristic of carotenoids having conjugated carbonyl groups 
[H, [HI- Addition of water into the stock solutions of 
astaxanthin in DMSO changes absorption spectra signif- 
icantly. All astaxanthin aggregates were prepared in a 
way that astaxanthin concentration, 20 /iM, remains con- 
stant in all solvent /water mixtures. 

In hydrated DMSO both J- and H-types of astaxan- 
thin aggregates are readily formed. At DMSO:water 
ratio of 1:1, a clear band 17450 cm -1 (573 nm) indi- 
cates that most of the astaxanthin molecules in solu- 
tion formed J-aggregates. At increasing water concen- 
tration, J-aggregates transforms to H-aggregates as in- 
dicated by narrow absorption band at 25640 cm -1 (390 
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nm). Thus, DMSO promotes formation of both types of 
aggregates with J-type generated only in rather narrow 
water concentration in agreement with earlier reports on 
astaxanthin [l8[ or zeaxanthin [l4[ in methanol. How- 
ever, while in previous studies using methanol or ace- 
tone required higher initial concentrations of carotenoid 
to induce stable J-aggregates, in DMSO the J-aggregates 
are generated even at moderate concentration of 20 fiM. 
Since we kept the concentration of astaxanthin constant 
during experiments, absorption spectra shown in Fig. Q] 
also demonstrate how the extinction coefficient of the 
S0-S2 transition changes upon aggregation. In both sol- 
vents, formation of aggregates is accompanied by signifi- 
cant decrease of extinction coefficient of both J- and H- 
aggregates. 



III. CALCULATION OF ABSORPTION 
SPECTRUM 

In order to characterize absorption spectra of astaxan- 
thin aggregates, we apply a line shape model based on 
second cumulant expansion of the electron-phonon cou- 
pling. In contrast to recently reported vibronic exciton 
approach aiming to explain absorption and CD spectra 
of aggregates of another carotenoid, lutein |l9j, we do not 
treat vibrational modes of the carotenoid explicitly. In 
the case of a carotenoid monomer, the two approaches 
would lead to the same results (second cumulant ex- 
pansion is exact for calculation of absorption spectra of 
harmonic systems The reduction to a purely elec- 

tronic exciton approach, which is very successful in de- 
scribing chlorophyll based aggregates [37], simplifies the 
treatment of the system significantly. It is also moti- 
vated by the lack of detailed vibrational structure in as- 
taxanthin and by only a small perturbation of the vi- 
brational band structure upon aggregation reported for 
other carotenoids [38] 

The astaxanthin molecule in solution can be repre- 
sented by the following Hamiltonian 



(e g +T + V g )\g)(g\ + (e e +T + V e )\e)(e\, (1) 



where the states \g) and |e) were introduced in Section 
HI e g and e e are the ground and excited state electronic 
energies, T represents the kinetic energy operator of the 
intramolecular nuclear modes of the carotenoid and the 
kinetic energy of the solvent DOF, and V g and V e rep- 
resent the potential energy surfaces of these DOF in the 
electronic ground and excited states, respectively. The 
Hamiltonian, Eq. can be split into the usual system 
(S), bath (B) and interaction (I) parts as follows 

ffcar = T+V^ + e g \g)(g\ + [e e + {V e -V g ) eq ]\e)(e[ 

Hb Hs 



(2) 




In the interaction Hamiltonian Hj we introduced so- 
called energy gap operator 



AV =V e -Vg- (V e - Vg) c 



(3) 



where (V e — V g ) oq — tr{(V e — V^)We q } is the equilibrium 
expectation value of the difference between the excited 
state and the ground state nuclear PES. 

To describe spectroscopy, let us introduce the light- 
matter interaction Hamiltonian in the semiclassical 
dipole approximation, i.e. 



H E {t) = -fi-E(t), 
where the dipole moment operator reads 
» = d\g)(e\+d*\e)(g\, 



(4) 



(5) 



with d independent of the solvent and intramolecular 
carotenoid DOF (Condon approximation). The absorp- 
tion coefficient can then be expressed in terms of linear 
response function S^'(t) as (33 ] 



a(w) 



n(w) 



Re / dco S (1) (t)e l 



(G) 



For a two-level system interacting with harmonic oscilla- 
tors, exact expression for the linear response function can 
written down in terms of the so-called line shape function 



l T 

9(t) = ^J dT J drC(r r ), 



(7) 



where C(t) is so-called bath or energy gap correlation 
function 



C(t) = tr{e* HBt AVe-* HBt AV W eq } 



(8) 



Here, W cq is the equilibrium density matrix of the in- 
tramolecular and solvent DOF. The first order response 
now reads as 



S<U(t) = !-Q{t)\d\ 2 e-^)- 



(9) 



The prefactor | results from the averaging of the response 
function over an isotropic distribution of the dipole mo- 
ment vectors of the molecule in space. By prescribing 
the correlation function C(t) in a form of two indepen- 
dent contributions 

C(t) = G,i b (t) + C solvent (t), (10) 

where the contribution of intramolecular vibrations 

Cvib(i) = A vib ( — — cos(w v ibi) 

- iu; vib sin(w vib t))e- r " b * /2 , (11) 
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Figure 2: Experimental (red) and calculated (black) absorp- 
tion spectrum of a monomer astaxanthin in DMSO. The 
parameters of calculation are: transition frequency u eg — 
20700 cm -1 , vibrational frequencies lj^I = 1520 cm -1 , and 
= 1150 cm -1 , reorganization energies of the vibrations 
A^ = 1220 cm -1 , and A^ = 250 cm -1 , damping time of the 
vibrations t^ 2 ' = = 1 ps, solvent reorganization energy 
A S oiv = 550 cm -1 , solvent correlation time r c = 100 fs. 

and the contribution of solvent reads as 

C solven t(*) = A solv (pj?- - iT^j e- Tt , (12) 

we completely specify the absorption spectrum of the 
monomer. Details of the model correlation functions 
Eqs. dTU and dHJ can be found in the Ref. (39j . 

The correlation function, Eq. (|5J), enables us to de- 
scribe another important contribution to the absorp- 
tion spectrum, so-called static disorder or inhomogeneous 
broadening. Provided that there is a certain distribu- 
tion of transition energies huj eg such that it is charac- 
terized by a constant energy gap correlation function 
Cdis.(t) = A 2 , this could be easily incorporated into the 
absorption spectrum formula, Eq. (jB]), which now reads 

oo 

a[uj) = i_I^r m / du e -*(t)-£f (13) 
3n(w) J 
o 

Instead of Eq. (flU)) , one can alternatively use explicit 
drawing of the energy gap value out of a Gaussian dis- 
tribution with the width A 2 . For a monomer, both ap- 
proaches are equivalent. For aggregate, the line-shape 
function depends on the aggregate structure, and the 
simulation of disorder has to be preformed explicitly. 

IV. EXCITONIC MODEL 

When aggregates are formed in the solution, astaxan- 
thin molecules start to interact with each other. We will 



use the subscript n and the superscript (n) to distinguish 
individual molecules of the aggregate in this section. The 
total Hamiltonian reads 

N 

H = Y,H&}+H int , (14) 

71=1 

where _ff; n t represents all possible interaction terms that 
follow from quantum chemistry. Electronic states of the 
aggregates can be expressed using the basis composed of 
excitations of individual monomers. The aggregate is in 
the ground state when all the molecules in the aggregate 
are in the ground state. Thus we can write 

N 

\G) = II \9n), (15) 

n=l 

where \g n ) is the ground state of the n th astaxanthin. 
The excited states of the aggregate in the spectral region 
around the 52 band of the monomeric astaxanthin are 
formed from single excitation states 

fc-l N 

\E k ) = J] \9n)\e k ) J] \9m), (16) 

n—1 m=k+l 

where \ek) is the excited S2 state of the k th astaxanthin 
in the aggregate. 

Expressed in these states, the aggregate Hamiltonian 
reads 

H = H B + ]T(e„ + AV^)\E n )(E n \ 

n 

+ Jmn\E m )(E n \, (17) 

where the ground state energy of the aggregate e g is set 

to zero, e„ = hJig is the excited state energy of the 11 th 

astaxanthin (e„ = ee+trbath{(V" e (n) -V s ( ™ ) )M^ e9 }) and J nm 
is the electrostatic coupling between two astaxanthins 

Jnm = trbath{(E n \Hi nt \E m )W eq }. (18) 

The electronic eigenstates of the Hamiltonian, 
Eq. (fTT|) , can be found by the diagonalization of its purely 
electronic part (i.e. ignoring AV( n) and H B in Eq. (fT7f ). 
The new excited eigenstates and their energies will be de- 
noted \E n ) and e n , respectively. The total Hamiltonian 
now reads 

H = H B + J^izn + AV nn )\E n )(E n \ 

n 

+ J2 &V mn \E m )(E n \. (19) 

The last term represents energy transfer between the elec- 
tronic eigenstates of the aggregate, the effects of which 
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lead to broadening of the absorption spectra. We will ig- 
nore these effects in further discussion and assume them 
to be smaller than the effects of broadening due to pure 
dephasing (the diagonal term AV nn ) and the disorder. 
The pure dephasing terms read as 

AV nn =Y,\(En\E k )\ 2 AV^. (20) 

k 

The delocalized eigenstates \E n ) are usually termed 
molecular excitons [40(. Thus we arrive at a situation 
which is similar to calculation of the absorption spec- 
trum of the monomer, Eq. (|13p . except of the number of 
involved states which is equal to the number of members 
of the aggregate. We thus have 

oo 

a(«) = ^^(V K| 2 Im / doj e-^~^~^) dlss , 
3n(u) ^ J 

(21) 

where (. . . }diss denotes explicit averaging over energy gap 
disorder, d n are the transition dipole moments of the n th 
eigenstate, and g nn {t) is the lineshape function corre- 
sponding to the energy gap operator AV nn in the Hamil- 
tonian Eq. (fl"9f . Assuming that the bath fluctuations on 
different molecules are the same, but are not correlated, 
we can write for the line shape functions 

9„n(t)=Y,\( E n\Ek)\ 4 g(t), (22) 
k 

where g(t) is the lineshape function of a monomeric 
molecule. 



partial charges. For preparation of the initial configura- 
tion, two astaxanthin molecules with distance of larger 
than 10 A were solvated in water box, which was fol- 
lowed by energy minimization for a few thousands steps. 
The system was then equilibrated for 2 ns followed by 
a 5 ns production run. The simulation was carried out 
in the NVT ensemble at 300 K, and temperature was 
controlled by the Berendsen thermostat [4g|. Equations 
of motion were integrated using the Leapfrog algorithm 
with the time step of 1 fs. The non-bonding interac- 
tions were cut at distance of 12.0 A, and long-range elec- 
trostatic interactions were treated by the Particle Mesh 
Ewald procedure [471 |48| . All bonds involving hydro- 
gen, atoms were constrained using the SHAKE algorithm 
[49j and a trajectory was constructed by sampling the 
coordinates at each 5 ps for further analysis such as dis- 
tance between conjugated chains and head groups. All 
MD simulations were performed by employing Amber85 
program package, and VMD program was used for vi- 
sualizations and preparation of snapshots and distances 
[50]. Fig. [3] shows typical snapshot from MD simulation 
which shows the aggregation taking place when astax- 
anthin molecules solvate in the aqueous solutions. The 
structure of astaxanthin dimer is stabilized within first 
few ps, which was revealed by analyzing the trajectory 
from MD simulation. The conjugated chains of two as- 
taxanthin molecules have in average the distance of 4.1 
A though the head groups can have even closer distances 
during the MD simulation. 



VI. RESONANCE COUPLING BETWEEN 
CAROTENOID MOLECULES 



V. MOLECULAR DYNAMICS SIMULATION 

In order to obtain parameters for calculation of ex- 
citonic interaction, aggregation of astaxanthin dimer in 
aqueous solution was studied by means of classical molec- 
ular dynamics (MD) simulations with non-polarizable 
force fields. The simulated system consisted of two as- 
taxanthin molecules and 2700 water molecules forming a 
slab of liquid in the center of the simulation box, with two 
water/vapor interfaces perpendicular to the z axis. The 
x, y and z-dimensions of the rectangular simulation box 
were set to 50.0, 50.0 and 100.0 A, respectively with pe- 
riodic boundary conditions in all three dimensions. The 
astaxanthin carotenoid molecule was modeled by using 
the general Amber force field (GAFF) [4l| and water 
molecules were modeled by using SPC/E model [42j. To 
obtain partial charges of astaxanthin molecule, ab initio 
geometry optimization of single astaxanthin molecule was 
performed using the Gaussian 03 package [43] by employ- 
ing the density functional theory B3LYP/cc-pVT. After 
the optimization by ab initio calculations, restricted elec- 
trostatic potential (RESP) procedure was applied using 
Antechamber program |44|,_which is implemented in the 
Amber program package [451 ] for obtaining the atomic 



Within the usual excitonic model one could calcu- 
late resonance coupling between two molecules in dipole- 
dipole approximation. Transition dipole moment vec- 
tor and mutual positions of the molecules are param- 
eters of the model. This approximation is applicable 
to well spaced compact molecules, such as chlorophylls 
in photosynthetic antennae (i.e. in proteins), but fails 
for chain-like molecules such as astaxanthins in solution. 
Obviously, the distance between two astaxanthins can be 
much smaller than the extent of the molecule, in which 
case dipole-dipole approximation breaks down. In order 
to estimate the coupling between astaxanthins, we need 
to use quantum chemical packages, or, as we do here, to 
construct a phenomenological model carotenoid with a 
qualitative description of the molecular states. We aim 
at a qualitative theory of coupling, which allows us to 
make claims beyond the validity of the dipole-dipole ap- 
proximation. 

The S*2 state of a carotenoid originates from the excita- 
tion of the conjugated 7r-electronic system. Therefore, we 
consider only the delocalized 7r-electrons, and we calcu- 
late its wavefunction using an Hiickel type semi-empirical 
method. Our model carotenoid is made of carbon atoms 
with bond length 1.4 A and with positions chosen in such 
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Figure 3: Structure of astaxanthin dimer resulting from 
molecular dynamics simulations. Carbon atoms are shown 
in blue, hydrogens are white and oxygens red. 



way, that all bonding angles are 120°. The structure is 
depicted on Fig. [TJ We describe the 7r-electron p-orbital 
wavefunctions by the STO-3G basis of atomic orbitals 

Inspired by the extended Huckel method (EHM) [52], 
we define the Huckel Hamiltonian in the form 



\H 



Huckel I 



Huckel 



6^a+(S^-6^)P, (23) 



where = {4>n\<i>v) are overlap integrals of the atomic 
orbitals \<p^) and a, j3 are parameters of the method con- 
structed from valence state ionization potentials [HH . By 
solving eigenvalue problem in a form 



E 



H, 



Huckel 



(24) 



we obtain the expansion coefficients Cj of the molecular 
orbitals in the basis of the atomic orbitals. Eq. (|24| can 
be rewritten as 



(25) 



with a 1 = a — j3 and = Si — ft. We fit the parameter 
a' to reproduce values of HOMO-LUMO transition ener- 
gies of several polyenes of a different length reported in 
Ref. [53j|. Because we are interested in only in the en- 
ergy difference, the particular value of the parameter /3 is 



irrelevant, and we may work with just a' . Unlike in the 
original EHM. this value does not have a direct physical 
interpretation. The comparison with polyens yields 



a 



66.48 eV. 



(26) 



Using the obtained molecular orbitals, we calculate the 
HOMO-LUMO energy gap and the related transition 
dipole moment in a standard manner. We get 



ae homo- 



LUMO 



HOMO-LUMO 



20742 cm" 
= 28.6 D . 



(27) 



In order to calculate the resonance coupling between 
two astaxanthins, we evaluate the corresponding matrix 
element of the electrostatic interaction Hamiltonian [54] 

V a h,a'V = J dx\... dr N y df i . . . df N 

{\ - 1 \ ^ Zj _ \ - Zi 
^Vi-H ^ki-Rjl ^IRi-^I 



V(ri,---,rjv) 



(28) 



Here, N is the number of all electrons of given molecule, 
Zi are the nuclear charges in units of the elementary 
charge, and are coordinates of electrons of astax- 
anthin 1 and 2, R/ and R/ are coordinates of nuclei of 
astaxanthin 1 and 2 and ^ a (ri, . . . , r^r), ^(f i, . . . , r^r) 
are complete electronic wavefunctions of astaxanthin 1 
and astaxanthin 2 respectively. They are constructed 
as Slater determinants from atomic orbitals. Overlaps 
between atomic orbitals of different astaxanthins are ne- 
glected. Contributions of particular terms of Hamilto- 
nian, Eq. (|28|) . can be found by using the Slater-Condon 
rules for matrix elements |511. 



VII. ABSORPTION SPECTRA OF 
ASTAXANTHIN AGGREGATES 

Molecular dynamics simulations show that in water 
carotenoids tend to minimize the exposure of their hy- 
drophobic chain to the solvent. Thus a parallel dimer 
(see Figs. [3] and |4j is a possible arrangement. The dipole- 
dipole coupling formula predicts a positive value J of the 
resonance coupling. Excitonic model predicts in such a 
dimer existence of two exited states, which have a form 
of a symmetric and anti-symmetric linear combinations 
of the S*2 excited states of individual astaxanthins. For 
the parallel dimer arrangement, optical excitation of the 
energetically higher state is allowed, whereas the lower 
lying state is dark. Thus one expects a blue shift of the 
absorption spectrum upon aggregation. 
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However, in a mixture of DMSO and water, at the 
percentage of water of 50 — 60 we observe a red shift 
(see Fig. 1). Such a shift requires negative value of the 
resonance coupling. Dipole-dipole coupling formula re- 
sults in negative coupling for a "serial" dimer (see Fig. 0]) 
where the centers of the astaxanthins are shifted by a 
sufficient distance X. As we discussed above, dipole- 
dipole coupling is not a reliable approximation for as- 
taxanthins, and we therefore performed a calculation for 
different values of X = 0, . . . , 30 A and R = 4.1 A (the 
value predicted by MD simulations in Section[Vj with our 
semi-empirical method. The resulting values of coupling 
are presented in Fig. [5] and the comparison of the two 
predictions is made in the inset. Both results coincide 
for large distances R and X where dipole approximation 
holds. One can immediately notice that dipole-dipole 
formula grossly overestimates the positive coupling for 
X = 0, it predicts a rapid switch to negative values of J 
(at X — 2.9 A) and even for larger X it overestimates the 
magnitude of the negative coupling. Our Hiickel-type cal- 
culation, on the other hand, predicts moderate values of 
the resonance coupling (between —1500 and 3000 cm" 1 ), 
and the switch from positive to negative values occur- 
ing at larger X, roughly when X is comparable with the 
half length of the astaxanthin conjugated chain. Since 
our calculations of the resonance coupling are too simple 
to aim at quantitative results, we will use the parame- 
ter J for selected structures as a fitting parameter in the 
next section. We also constrain the resonance coupling 
to values up to few thousands of cm . 



4000 



a) Parallel dimer 



J> 



b) "Serial" dimer 




J < 



c) chain "serial" multimer 



d) parallel trimer 



J < 




J>0 



Figure 4: Characteristic structures of the astaxanthin dimers 
and aggregates used for calculation of absorption spectra. 
Part a) Parallel dimer and similar slightly shifted arrange- 
ments are characterized by J > 0. Part b) The "serial" dimer 
with large shift X is characterized by J < 0. Part c) Larger 
aggregates of the "serial" character can have a form of chains. 
Part d) The parallel dimer can form stacks or more compact 
structures such as a trimer with equal distance between the 
conjugated chains of the three molecules. 

To analyze the absorption spectra of astaxanthin in 
DMSO with varying percentage of water, we calculated 
spectra of the two representative mutual configurations 



— dipole-dipole 

— full-electrostatic 




-2000 



10 15 20 

Shift X [A] 



25 



30 



Figure 5: Dependence of the resonance coupling on the mu- 
tual shift X (see Fig. [4j of the astaxanthin molecules . 



of astaxanthin molecules, namely those with positive and 
negative resonance coupling energy J. We start with a 
simple idea that with increasing percentage of water in 
the solutions, it will become energetically advantageous 
for the molecules to minimize the surface of their hy- 
drophobic chains exposed to water. With small amount 
of water, the free energy profile will be rather flat and 
we might preferentially find the molecules in the arrange- 
ment similar to that one of the "serial" dimer (see Fig. 0} . 
With increasing amount of water a strict parallel dimer 
arrangement will be more and more preferred. 

It is in no way plausible that one could reproduce the 
absorption spectrum by composing it from a small num- 
ber discrete structures, e.g. dimers and trimers of partic- 
ular value of coupling J. It is more likely that there exists 
a smooth distribution of structures depending on the par- 
ticular interactions between molecules and the solvent at 
a given concentration of water. This would result in a dis- 
tribution of values of the coupling J. A direct simulation 
of the disordered structure of aggregates requires some 
prescription of the probability of various configuration. 
In the absence of a physical model supporting the distri- 
bution, one might attempt an effective approach, where 
the fluctuation of the structure would be simulated by 
fluctuations of the transition energies of the carotenoids 
(disorder parameter A) in several selected structures. If 
successful, it could provide guidance in constructing a 
more detailed model of aggregation which would take into 
account the thermodynamics of the process. 

The exciton model is not able to account for the re- 
duced extinction coefficient of the aggregates. The total 
calculated extinction coefficient of the aggregates would 
thus be overestimated. As a result, one cannot expect to 
deduce actual concentrations of the aggregates from our 
model, because the transition dipole moment of an aggre- 
gate is reduced by an unknown factor. In order to fit the 
experimental spectra, we will therefore calculate normal- 
ized absorption spectra a\{uj), 02 (w, J2), a(w, J3),. . . for 
monomer, dimer, trimer and larger aggregates with some 



values J of the resonance coupling. We will then combine 
them to get the total normalized absorption spectrum 
a{uS) as 

1 K 

ot{u) = jj ^2 a k a k(u, Jfc), (29) 

n— 1 

where J\f is a normalization constant. This will give us 
an qualitative picture of the ratio of various aggregates 
in the mixture, 

The values of J and the selected structures where cho- 
sen according to the features observed in the astaxanthin 
spectra at 50 % and 90 % of water (see Fig.rTJV) while tak- 
ing into account qualitative results of the molecular or- 
bital simulations. Such a procedure cannot lead to unique 
results, but still seriously constrains the composition of 
the aggregate spectra, as demonstrated below. 



VIII. DISCUSSION 

A. J-Aggregate Spectrum 

Let us first analyze the spectra measured at 50 % of 
water. The experimental spectrum shows a significant 
red shift by w 2800 cm -1 . This is consistent with for- 
mation of J-aggregates (J < 0) i.e. with the "serial" 
dimer and larger aggregates of the similar form. Combi- 
nation of the monomeric spectrum with a spectrum of 
the dimer (J = —2500 cm -1 ) leads to the spectrum 
with characteristic shoulder. Adding small portions of 
trimers, tetramers and larger aggregates up-to six mem- 
bers leads to a good fit of the experimental spectrum 
(black dashed curve in Fig. |6j) starting from 16000 cm , 
with slightly insufficient absorption for frequencies larger 
than 20000 cm -1 . The spectral region between 14000 
cm -1 and 16000 _1 probably contains aggregates with 
even more members than six. By adding the parallel 
dimer spectrum (taken from the fit of the spectra at 90 
% of water, see below) we obtain far better fit for the 
frequencies above 20000 cm -1 confirming the presence of 
some parallel H-aggregate type dimers already at 50 % 
of water. 

The parameters of the fit are summarized in Tab. fl] 
The fit is composed of absorption spectra of small aggre- 
gates whose constituents have the parameters same as 
the monomer in Fig. [2] . The influence of the trimer and 
larger aggregate spectra on the total spectrum is small 
and it was not possible to make conclusions on the mag- 
nitude of the disorder in them. We leave their disorder 
the same as in the monomer. 



B. H- Aggregate Spectrum 

At the concentration of water equal to 90 %, we can as- 
sume the formation of parallel dimers and possibly larger 




Figure 6: Absorption spectrum of a mixture of astaxanthin 
aggregates at 50 % of water. Theoretical spectrum is com- 
posed of aggregates up to six members with J < (parallel 
dimer) and of a sizable amount of parallel dimers with J > 
(see Tab. [1} . The spectra of various aggregates are all normal- 
ized for their maximum equal to 0.3. 



Nr. of members 


A [cm" 1 ] 


J [cur 1 ] 


w 


1 


800 




0.50 


2 


800 


-2500 


0.75 


3 


800 


-2500 


0.28 


4 


800 


-2500 


0.09 


5 


800 


-2500 


0.03 


6 


800 


-2500 


0.02 


2 


2400 


2100 


0.21 



Table I: Parameters of the fit in Fig. [6] First column repre- 
sents the number of monomers in an aggregate. The transition 
energies huj eg are those of the monomer in Fig. [2] (parameters 
in the text) with disorder parameter A and nearest neighbor 
coupling J. The sum spectrum is composed of the normal- 
ized spectra of various aggregates according to Eq. (|29[) with 
coefficients a n . The last row is the parallel dimer from Tab.HTl 



aggregates of similar structure. Positive resonance cou- 
pling is suggested by the blue shift of the spectrum with 
respect to the spectrum of the monomer. The line shape 
can be fitted by a mixture of monomers, dimers and 
trimers of the parallel type, and a small contribution of 
the "serial" dimers. The parallel dimer has to have a sig- 
nificantly larger disorder than the monomer, suggesting 
the existence of coupling disorder or some deformation of 
the chain . The trimer contribution cannot be explained 
by simple stacking of the monomers. This would result in 
a nearest neighbor coupling and a second nearest neigh- 
bor coupling significantly reduced, and correspondingly 
in an insufficient shift of the trimeric peak. We have 
to rather assume an increased excitonic coupling which 
acts among all three molecules of the aggregate. This is 
consistent with all three molecules touching by their con- 
jugated chains (see Fig.fJJ)) and a compactification of the 
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aggregate (molecules are closer to each other than in the 
dimer). As an alternative to this fit, one could try to ex- 
plain the most blue peak of the experimental H-aggregate 
spectrum of astaxanthin (Fig. [7]) by a dimer with large 
resonance coupling. This coupling would, however, have 
to be almost an order of magnitude larger than the one 
calculated by our semi-empirical method. Assuming a 
compactified trimer returns the required values of reso- 
nance coupling into a reasonable parameter region. 



In both investigated absorption spectra (Figs. [5] and 
[7]) we can see several distinct excitonic features. The 
most significant feature is the spectral shift of the bands. 
However, the line shapes also bear signatures of band nar- 
rowing, evident e.g in the appearance of the vibrational 
shoulder in the dimer spectrum at 50 % of water, or in 
the narrow line shape of the trimers at 90 % of water. 
Both these shapes determine the most obvious features 
of the aggregate spectra. 




Figure 7: Absorption spectrum of a mixture of astaxanthin 
aggregates at 90 % of water. The spectrum is composed of 
mono-, di- and trimers with J > (parallel dimer) and of 
a small amount of "serial" dimers with J < (see Tab. ITTj) . 
The spectra of various aggregates are all normalized for their 
maximum equal to 0.3. 



Nr. of members 


A [cm" 1 ] 


J [cm" 1 ] 


w 


1 


800 




0.36 


2 


2400 


2100 


0.87 


3 


1200 


2750 


0.60 


2 


800 


-2500 


0.06 



Table II: Parameters of the fit in Fig. [Tl] First column repre- 
sents the number of monomers in an aggregate. The transition 
energies fkd eg are those of the monomer in Fig. [2] (parameters 
in the text) with disorder parameter A. In case of a dimer, J 
represents the coupling between the two monomers, for trimer 
J is the coupling of a monomer to other two monomers. The 
sum spectrum is composed of the normalized spectra of par- 
allel mono-, di- and trimers and "serial" dimer. The sum 
spectrum is composed of the normalized aggregate spectra 
according to Eq. (|29[) with coefficients a n . The last row is the 
"serial" dimer from Tab. [T] 



IX. CONCLUSIONS 

We have measured and analyzed absorption spectra 
of aggregates of the carotenoid astaxanthin in hydrated 
dimethylsulf oxide. We have demonstrated that these 
absorption spectra can be explained quantitatively as 
a sum of spectra of small aggregates. Phenomenolog- 
ical excitonic model of resonance interaction between 
7T— conjugated chains of astaxanthin molecules was con- 
strained by semi-empirical calculations. The observed 
spectral shifts were assigned to formation of J-type and 
H-type aggregates and their changing ratio under differ- 
ent concentration of water. Fitting experimental absorp- 
tion spectra lead to a qualitative estimation of ratio of 
aggregates of various types in the sample at 50 % and 
90 % of water in dimethylsolfoxide. With increased per- 
centage of water in the solvent, the proportion of compact 
aggregates with a decreasing surface area exposed to the 
solvent is growing. Based on absorption spectra, it is 
not possible to exclude that different aggregate configu- 
rations result in the same observed features. The ease 
with which the spectra can be explained by the spectra 
of small aggregates, however, provides a relatively high 
degree of confidence in our results. 



Acknowledgments 

The research was supported by grants from the Czech 
Ministry of Education (MSM6007665808) and the Czech 
Science Foundation (202/09/1330). J.O. acknowledges 
the support of GAUK through grant nr. 416311. The au- 
thors thank Lucie Tesnohlidkova for her help with prepa- 
ration of astaxanthin aggregates. 



[1] T. Polivka and V. Sundstrom, Chem. Rev. 104, 2021 
(2004). 

[2] E. Ostroumov, M. Muller, C. Marian, M. Kleinschmidt, 
and A. Holzwarth, Phys. Rev. Lett. 103, 108302 (2009). 



[3] M. Maiuri, D. Polli, D. Brida, L. Luer, A. LaFountain, 
M. Fuciman, R. J. Cogdell, H. Frank, and G. Cerullo, 
Phys. Chem. Chem. Phys. 14, 6312 (2012). 

[4] T. Polivka and V. Sundstrom, Chem. Phys. Lett. 477, 1 



10 



(2009) . 

[5] D. Kosumi, T. Kusumoto, R. Fujii, M. Sugisaki, Y. 
Iinuma, N. Oka, Y. Takaesu, T. Taira, M. Iha, H. A. 
Frank, and H. Hashimoto, Phys. Chem. Chem. Phys. 22, 
10762 (2011). 

[6] C. C. Gradinaru, J. Kennis, E. Papagiannakis, I. van 
Stokkum, R. Cogdell, G. Fleming, R. Niederman, and R. 
van Grondelle, Proc. Natl. Acad. Sci. U.S.A. 98, 2364 

(2001) . 

[7] T. Polivka and H. A. Frank, Acc. Chem. Res. 43, 1125 

(2010) . 

[8] A. Ruban, R. Berera, C. Ilioaia, van Stokkum I.H.M., 

J. Kennis, A. Pascal, H. van Amerongen, B. Robert, P. 

Horton, and R. van Grondelle, Nature 450, 575 (2007). 
[9] Y. Koyama, F. S. Rondonuwu, R. Fujii, and Y. Watan- 

abe, Biopolymers 74, 2 (2004). 
[10] J. H. Starcke, M. Wormit, J. Schirmer, and A. Dreuw, 

Chem. Phys. 329, 39 (2006). 
[11] D. Ghosh, J. Hachmann, T. Yanai, and G. K. L. Chan, 

J. Chem. Phys. 128, 144117 (2008). 
[12] M. Kleinschmidt, M. M. C, M. Waletzke, and S. Grimme, 

J. Chem. Phys. 130, 044708 (2009). 
[13] M. Simonyi, Z. Bikadi, F. Zsila, and J. Deli., Chirality 

15, 680 (2003). 

[14] H. H. Billsten, V. Sundstrom, and T. Polivka, J. Phys. 

Chem. A 109, 1521 (2005). 
[15] C. Kopsel, H. Moltgen, H. Schuch, H. Auweter, K. Klein- 

ermanns, H.-D. Martin, and H. Bettermann, J. Mol. 

Struct. 750, 109 (2005). 
[16] Y. Mori, K. Yamano, and H. Hashimoto, Chem. Phys. 

Lett. 254, 84 (1996). 
[17] A. V. Ruban, P. Horton, and A. J. Young, J. Photochem. 

Photobiol. B 21, 229 (1993). 
[18] R. Giovannetti, Alibabaei, L. Pucciarelli, and F., Spec- 

trochim. Acta A 73, 157 (2009). 
[19] F. Spano, J. Am. Chem. Soc. 131, 4267 (2009). 
[20] C. Wang and T. M.J., J. Am. Chem. Soc. 132, 13988 

(2010). 

[21] C. Wang, D. E. Schlamadinger, V. Desai, and M. J. 
Tauber, Chem. Phys. Chem. 12, 2891 (2011). 

[22] A. Sujak, P. Mazurek, and W. I. Gruszecki, J Photochem. 
Photobiol. B 68, 39 (2002). 

[23] W. Gruszecki, in Photochemistry of Carotenoids, edited 
by G. B. H. A. Frank, A. J. Young and R. J. Cogdell 
(Kluwer Academic Publishers, Dordrecht, Netherlands, 
1999), Chap. Carotenoids in membranes, p. 363. 

[24] M. Aspinall-O'Dea, M. Wentworth, A. Pascal, B. Robert, 
A. Ruban, and P. Horton, Proc. Natl. Acad. Sci. U.S.A. 
99, 16331 (2002). 

[25] P. Chabera, M. Durchan, P. M. Shih, C. A. Kerfeld, and 
T. Polivka, BBA-Bioenergetics 1807, 30 (2011). 

[26] F. G. Gao, A. J. Bard, and L. D. Kispert, J. Photochem. 
Photobiol. A 130, 49 (2000). 

[27] J. Pan, G. Benko, Y. H. Xu, T. Pascher, L. Sun, V. Sund- 
strom, and T. Polivka, J. Am. Chem. Soc. 124, 13949 

(2002) . 

[28] J. Pan, Y. Xu, L. Sun, V. Sundstrom, and T. Polivka, J. 

Am. Chem. Soc. 126, 3066 (2004). 
[29] J. Xiang, F. S. Rondonuwu, Y. Kakitani, R. Fujii, Y. 

Watanabe, Y. Koyama, H. Nagae, Y. Yamano, and M. 

Ito, J. Phys. Chem. B 109, 17066 (2005). 
[30] X. F. Wang, Y. Koyama, H. Nagae, Y. Yamano, M. Ito, 

and Y. Wada, Chem. Phys. Lett. 420, 309 (2006). 
[31] J. Smith, M.B.and Michl, Chem. Rev. 110, 6891 (2010). 



[32] M. Cianci, R. P.J., A. Olczak, J. Raftery, N. Chayen, P. 

Zagalsky, and J. Helliwell, Proc. Natl. Acad. Sci. U.S.A. 

15, 9795 (2002). 
[33] A. van Wijk, A. Spaans, N. Uzunbajakava, C. Otto, H. 

De Groot, J. Lugtenburg, and F. Buda, J. Am. Chem. 

Soc. 127, 1438 (2005). 
[34] F. H. A., J. A. Bautista, J. Josue, Z. Pendon, R. G. Hiller, 

F. P. Sharpies, D. Gosztola, and M. R. Wasielewski, J 
Phys Chem B 104, 4569 (2000). 

[35] D. Zigmantas, R. G. Hiller, F. P. Sharpies, H. A. Frank, 
V. Sundstrom, and T. Polivka, Phys Chem Chem Phys 
6, 3009 (2004). 

[36] R. Doll, D. Zueco, M. Wubs, S. Kohler, and P. Hanggi, 

Chem. Phys. 347, 243 (2008). 
[37] M. H. Cho, H. M. Vaswani, T. Brixner, J. Stenger, and 

G. R. Fleming, J. Phys. Chem. B 109, 10542 (2005). 
[38] C. Wang, C. J. Berg, C.-C. Hsu, B. A. Merril, and M. J. 

Tauber, J. Phys. Chem. B article ASAP (2012). 

[39] S. Mukamel, Principles of nonlinear spectroscopy (Oxford 
University Press, Oxford, 1995). 

[40] H. van Amerongen, L. Valkunas, and R. van Gron- 
delle, Photo synthetic Excitons (World Scientific, Singa- 
pore, 2000). 

[41] J. Wang, R. M. Wolf, J. W. Caldwell, P. A. Kollman, and 
D. A. Case, J. Comput. Chem. 25, 1157 (2004). 

[42] H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, 
J. Phys. Chem. 91, 6269 (1987). 

[43] M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuse- 
ria, M. A. Robb, J. R. Cheeseman, J. Montgomery, J. A., 
T. Vreven, K. N. Kudin, J. C. Burant, J. M. Millam, S. S. 
Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, 
G. Scalmani, N. Rega, G. A. Petersson, H. Nakatsuji, M. 
Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, 
M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, 
M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J. B. 
Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, 
R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, 
C. Pomelli, J. W. Ochterski, P. Y. Ayala, K. Morokuma, 

G. A. Voth, P. Salvador, J. J. Dannenberg, V. G. Za- 
krzewski, S. Dapprich, A. D. Daniels, M. C. Strain, O. 
Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, 
J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clif- 
ford, J. Cioslowski, B. B. Stefanov, G. Liu, A. Liashenko, 
P. Piskorz, I. Komaromi, R. L. Martin, D. J. Fox, T. 
Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, 
M. Challacombe, P. M. W. Gill, B. Johnson, W. Chen, 
M. W. Wong, C. Gonzalez, and J. A. Pople, Gaussian, 
Inc., Wallingford CT, 2004. 

[44] J. Wang, W. Wang, P. A. Kollman, and D. A. Case, J. 

Mol. Graph. Model 25, 247260. (2006). 
[45] D. A. Case, T. A. Darden, T. E. I. Cheatham, C. L. Sim- 

merling, J. Wang, R. E. Duke, R. Luo, K. M. Merz, B. 

Wang, D. A. Pearlman, M. Crowley, S. Brozell, V. Tsui, 

H. Gohlke, J. Mongan, V. Hornak, G. Cui, P. Beroza, C. 
Schafmeister, J. W. Caldwell, W. R. Ross, and P. A. Koll- 
man, Amber 8, University of California, San Francisco, 
2004. 

[46] H. J. C. Berendsen, J. P. M. Postma, W. F. van Gun- 
steren, A. DiNola, and J. R. Haak, J. Chem. Phys. 81, 
3684 (1984). 

[47] T. Darden, D. York, and L. G. Pedersen, J. Chem. Phys. 

98, 10089 (1993). 
[48] U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, 

H. Lee, and L. G. Pedersen, J. Chem. Phys. 103, 8577 



11 



(1995). 

[49] J. P. Ryckaert, G. Ciccotti, and H. J. C. Berendsen, J. 
Comput. Phys. 23, 327 (1977). 

[50] W. Humphrey, A. Dalke, and K. Schulten, J. Mol. Graph- 
ics 14, 33 (1996). 

[51] A. Szabo and N. S. Ostlund, Modern Quantum Chem- 
istry Introduction to Advanced Electronic Structure The- 
ory (Dover, New York, 1996). 



[52] R. Hoffmann, J. Chem. Phys. 39, 1397 (1963). 

[53] S. M. Smith, A. N. Markevitch, D. A. Romanov, X. Li, 

R. J. Levis, and H. B. Schlegel, J. Phys. Chem. A 108, 

11063 (2004). 

[54] M. E. Madjet, A. Abdurahman, and T. Renger, J. Phys. 
Chem. B 110, 17268 (2006). 



